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EXECUTIVE  SUMMARY 


OBJECTIVE 

The  objective  is  to  document  and  analyze  two  weeks  of  surface  and  upper-air  meteorological 
measurements  made  in  support  of  an  international  effort,  coordinated  by  the  NATO/AC23/Panel  3 
Research  Study  Groups  8  and  21,  to  assess  the  effects  of  the  lower  atmosphere  on  the  propagation 
of  an  electromagnetic  phase  front  at  frequencies  from  10  to  94  GHz. 

RESULTS 

Surface  and  upper-air  meteorological  measurements  were  made  at  Gavres,  France  (47.42  N, 
3.2  W).  A  total  of  20  radiosondes  were  launched  on  seven  different  days.  In  all  cases,  the  dominant 
propagation  phenomenon  was  the  evaporation  duct. 

RECOMMENDATION 

It  is  strongly  recommended  that  phase-front  angle-of-arrival  estimates  calculated  from  surface 
meteorological  conditions  be  compared  to  measured  values. 

CONCLUSIONS 

The  numerical  modeling  of  propagation  from  waveguide  theory  at  10.5  and  35  GHz  in  normal 
and  ducting  environments  indicates  that  the  evaporation  duct  can  raise  the  angle-of-arrival  from 
about  1.5  mrad  to  about  6.2  mrad.  In  addition,  surface  roughness,  at  least  on  short  propagation  paths, 
has  little  or  no  effect  on  angle-of-arrival. 
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INTRODUCTION 

The  propagation  of  electromagnetic  (EM)  waves  is  critically  dependent  on  atmospheric  condi¬ 
tions.  At  millimeter  wavelengths,  one  of  the  most  well-known  propagation  effects  is  molecular 
absorption.  It  is  caused  by  the  interaction  of  the  passing  EM  wave  and  the  charge  distribution  of  a 
molecule.  In  the  troposphere,  the  major  effects  of  molecular  absorption  are  attributable  to  the  gas¬ 
eous  constituents  of  O2  and  H2O  molecules.  Although  molecular  absorption  is  generally  thought  to 
be  the  dominant  propagation  mechanism,  Anderson  [1990]  demonstrated  that  the  evaporation  duct, 
a  persistent  propagation  phenomenon  over  water,  is  capable  of  trapping  the  EM  energy  and  guiding 
it  along  the  earth’s  surface  to  ranges  of  more  than  twice  the  radio  horizon  as  if  there  were  no  molecu¬ 
lar  absorption.  These  measurements,  made  at  94  GHz  on  a  40-km  overwater  path,  showed  that  the 
evaporation  duct  typically  enhanced  the  received  signal  strength  by  some  60  dB,  which  almost  com¬ 
pensated  for  the  loss  due  absorption. 

At  short  ranges  and  low  terminal  heights,  well  within  the  radio  horizon,  the  evaporation  duct  has 
been  shown  to  be  a  significant  propagation  phenomenon  at  X-band  (Dockery,  1988;  Anderson, 
1993).  Anderson’s  radar  measurements  {op.  cit.)  clearly  indicate  that  the  vertical  modified  refractive 
index  profile  describing  the  evaporation  duct  must  extend  from  the  surface  through  the  mixed  layer, 
essentially  encompassing  the  entire  boundary  layer.  The  refractive  gradient  above  the  surface  layer 
has  a  substantial  effect  on  propagation  at  moderate  ranges.  Typically  it  affects  propagation  from 
about  50%  to  100%  of  the  horizon  range. 

The  low-altitude  propagation  of  millimeter  wavelength  signals  over  short  ranges  is  expected  to 
be  dominated  by  both  absorption  and  evaporation  ducting.  Surface  roughness  effects  (Miller, 
Brown,  and  Vegh,  1984)  may  also  be  a  significant  factor  in  the  propagation  of  millimeter  waves. 
For  targeting  and  tracking  systems  that  make  use  of  the  angle-of-arrival  of  a  ray  (direction  of  the 
normal  to  a  wave  front),  ducting  and  surface  roughness  effects  are  of  primary  concern.  Absorption 
reduces  the  signal  strength  but  does  not  change  the  dire^  tion  of  a  ray  as  ducting  and  scattering  do. 

The  objective  of  a  measurement  and  analytical  effort  sponsored  by  NATO/AC243/Panel  3 
Research  Study  Groups  8  and  21  was  to  characterize  surface  and  atmospheric  effects  on  the  propaga¬ 
tion  of  EM  phase  fronts  at  millimeter  wavelengths.  The  participants  are  listed  in  table  1.  Dutch, 
French,  and  German  groups  provided  EM  transmitter  and  receiver  systems  using  array  or  monopulse 
antennas  to  measure  the  angle-of-arrival.  The  Dutch  systems  (FEL/TNO)  operated  at  10.5  and  16 
GHz;  the  French  systems  (CELAR  and  CERT)  operated  at  35  and  36  GHz;  the  German  systems  oper¬ 
ated  at  35  and  94  GHz.  All  five  participants  made  surface  meteorological  measurements  (CELAR 
and  FEL/TNO  had  buoys  midway  between  He  de  Groix  and  Gavres).  DREV  and  NRaD  with 
NAVEURMETOCCEN  made  upper-air  observations  using  radiosondes. 

The  participants  agreed  to  document  their  particular  measurements  and  share  data  on  a  case-by¬ 
case  basis.  In  fulfillment  of  this  agreement,  the  measurement  program  will  be  reviewed  briefly  in 
this  report  (greater  detail  will  be  available  from  the  other  national  reports).  First,  the  physics  of  the 
surface  layer  meteorological  models  will  be  examined;  next,  surface  and  upper-air  measurements 
made  by  NRaD  and  NAVEURMETOCCEN  will  be  presented;  and,  finally,  some  calculations  of 
expected  ray  angles-of-arrival  will  be  analyzed. 
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Table  1 .  List  of  participants. 


Country 

Participant 

Canada 

Defense  Research  Establishment  Valcartier  (DREV) 

France 

Centre  d’Electronique  de  TArmement  (CELAR)  Centre  d’Etudes  et  de 
Recherches  de  Toulouse  (CERT) 

Germany 

Forschungsgesellschaft  fuer  Angewandte  Naturwissenschanften- 
Forschungsinstitut  fur  Hochfrequenzphysik  (FGAN-FHP) 

Netherlands 

Fysisch  en  Elektronish  Lab.  TNO  (FEL/TNO) 

United  States 

Naval  Command,  Control  arid  Ocean  Surveillance  Center, 

RDT&E  Division  (NRaD) 

Naval  European  Meteorology  and  Oceanography  Center 
(NAVEURMETOCCEN) 

MEASUREMENT  PROGRAM 


OVERVIEW 


Figure  1  shows  the  geographical  region  near  Lorient,  France,  in  the  Bay  of  Biscay.  A  9.8-km 
overwater  propagation  path  extends  from  He  de  Groix  to  Gavres  and  is  shown  by  the  dashed  line. 
The  transmitters  were  mounted  approximately  14.8  meters  above  mean  sea  level  on  He  de  Groix  and 
the  receivers  were  mounted  between  6  and  10  meters  above  mean  sea  level  at  a  site  near  Gravres. 
Tidal  variations  in  this  area  are  about  2  meters.  At  high  tide,  the  radio  horizons  are  14.7  km  for  the 
transmitter  and  8.2  and  11.6  km,  respectively,  for  the  low-  and  high-sited  receiver  antennas.  Even 
at  high  tide,  the  propagation  path  is  well  within  the  radio  horizon.  Winds  blowing  from  the  east- 
southeast  (about  140  degrees)  through  due  west  (270  degrees)  are  considered  to  be  from  the  sea.  With 
winds  from  these  directions,  it  is  assumed  that  the  meteorological  measurements  made  at  the  receiv¬ 
er  site  (near  Gavres)  are  representative  of  conditions  along  the  propagation  path.  Winds  blowing 
from  northwest,  through  north,  through  east-southeast  are  considered  continental  winds.  These  con¬ 
tinental  winds  indicate  that  the  meteorological  measurements  at  the  receiver  site  are  probably  not 
representative  of  conditions  along  the  propagation  path. 


Table  2  lists  the  data  and  time  of  each  upper-air  sounding  that  was  made  by  NRaD  and  N  AVEUR- 
METOCCEN.  Surface  meteorological  observations  of  pressure,  air  temperature,  relative  humidity, 
wind  speed,  wind  direction,  and  sea  surface  temperature  are  listed  in  this  table.  These  observations 
were  made  typically  a  few  minutes  before  the  radiosonde  launch.  Two  columns,  one  labeled  “Sur¬ 
face  Roughness"  and  the  other  labeled  “Monin-Obukhov  Length,”  list  important  parameters  that 
are  derived  from  considerations  of  the  energy  balance  in  the  surface  layer.  These  parameters  are 
examined  in  the  following  section.  The  last  three  columns  in  this  table  list  estimates  of  the  modified 
refractivity  gradient  with  respect  to  height.  These  gradients  are  discussed  in  subsequent  sections. 
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NAVEURMETOCCEN  provided  a  two-man  Mobile  Environmental  Team  (MET)  with  equip¬ 
ment  and  sensors  to  make  both  surface  and  upper-air  meteorological  measurements.  The  upper-air 
equipment  consisted  of  a  ;>tandard  US  Navy  Mini  Rawin  System  (MRS)  and  Vaisala  RS-80  radio¬ 
sondes  with  a  wind  velocity  and  direction  capability.  The  MRS  is  a  sophisticated  receiver  and  pro¬ 
cessing  system.  It  receives  RF  signals,  at  403  MHz,  transmitted  by  the  radiosonde  as  a  helium-filled 
balloon  carries  the  radiosonde  through  the  atmosphere.  The  radiosonde  instrument  package  directly 
senses  the  pressure,  temperature,  and  relative  humidity  of  the  air  parcel  and  telemeters  this  informa¬ 
tion  uack  to  the  ground-based  MRS  receiver  as  raw  data  (a  pulse-width  modulated  signal).  In 
addition,  it  retransmits  a  NAVAID  signal  (typically  OMEGA)  back  to  the  MRS,  which  processes 
this  retransmitted  data  to  determine  wind  direction  and  speed.  Internally,  the  MRS  monitors  the  raw 
pressure,  temperature,  and  humidity  data  stream  (a  new  value  is  sent  about  once  every  1 .3  seconds), 
converts  the  raw  data  into  true  units  (mb,  °C,  and  percent).  It  also  computes  the  height  and  refractiv- 
ity.  Further,  the  processor  smooths  the  data,  removes  outliers,  and  provides  a  number  of  display  and 
ouq>ut  options. 

NRaD  upper-air  equipment  consisted  of  a  Vaisala  model  PP- 1 1  processing  system,  which  is  less 
sophisticated  than  the  MRS.  The  PP-1 1  receives  and  decodes  the  same  raw  pressure,  temperature, 
and  humidity  data  stream  as  the  MRS,  but  the  PP-1 1  simply  converts  the  data  (to  true  units  of  mb, 
°C.  and  percent)  and  outputs  the  converted  stream  as  a  set  of  ASCII  characters.  There  is  no  smooth¬ 
ing  or  any  additional  processing  of  the  data. 

Both  the  MRS  and  the  PP- 1 1  read  a  punched-pc^  tape  to  get  a  set  of  coefficients  that  are  used 
to  convert  the  raw  data  stream  into  true  units.  Each  radiosonde  is  calibrated  at  the  factory  and  has 
its  own  punched  tape.  However,  the  PP-11  and  MRS  conversion  algorithms  are  identical. 

Although  less  than  half  the  observations  may  be  considered  representative  of  conditions  along 
the  propagation  path,  there  is  one  case  where,  for  two  days  in  a  row,  the  winds  were  from  the  sea. 
Therefore,  the  analysis  will  concentrate  on  measurements  made  on  these  two  days — 20  and  2 1  Sep¬ 
tember.  However,  before  proceeding  into  the  data  analysis,  the  meteorological  models  will  be 
quickly  reviewed  (greater  details  are  readily  available  in  the  open-literature  citations). 

METEOROLOGICAL  MODEL  OF  THE  SURFACE  LAYER 

Radio  wave  propagation  in  the  troposphere  is  controlled  by  the  spatial  distribution  of  the  radio 
refractive  index,  n,  and  by  the  surface  characteristics.  A  numerically  convenient  term  is  the  refractiv- 
ity  N  that  is  related  to  the  refractive  index  by  N  =  (n  -  1)10®.  The  refractivity  of  a  parcel  of  air  is 

n  =  +  (1) 

where  T  is  the  temperature  in  Kelvin,  P  is  the  pressure  in  mb,  and  E  is  the  water  vapor  pressure  in 
mb.  For  a  standard  or  normal  atmosphere,  the  refr{u:tivity  monotonically  decreases  with  increasing 
height  at  a  rate  of  39  ATkm  [Bean  and  Dutton,  1966].  By  Snell’s  Law.  a  ray  (direction  of  a  normal 
to  the  wave  front)  initially  launched  at  an  angle  parallel  to  the  earth’s  surface  will  be  refracted  down¬ 
ward  as  it  travels  in  range.  The  curvature  of  the  ray  path  is  less  than  the  curvature  of  the  earth,  so, 
as  the  ray  travels  in  range,  the  height  of  the  ray  above  the  surface  increases. 

The  modified  refractive  index,  M,  is  a  mapping  of  the  vertical  refractivity  from  a  curved  geome¬ 
try  to  Cartesian  coordinates  and  greatly  reduces  the  computational  complexity  of  radio  propagation 


4 


TD  1668 


analysis.  It  is  defined  as  M  =  N  +  {Z/a)10^,  where  Z  is  the  height  above  the  earth  's  surface  and  a  is 
the  earth’s  radius  (6371  km). 

In  the  troposphere,  measurements  of  T,  P,  and  RH  are  sufficient  to  compute  the  refractivity  of 
the  air  parcel.  However,  close  to  the  surface,  the  meteorological  measurements  are  strongly  in¬ 
fluenced  by  turbulence.  Instead  of  using  the  instantaneous  measurements  of  T,  P,  and  RH  to  charac¬ 
terize  the  refractivity  close  to  the  surface  (in  the  surface  layer),  the  accepted  approach  is  to  model 
the  turbulent  transport  mechanisms  and  infer  the  refractivity. 

While  there  is  no  exact  definition  of  surface  layer,  qualitatively,  it  is  that  part  of  the  atmosphere 
immediately  above  the  surface  where  the  fluxes  of  momentum,  heat,  and  moisture  can  be  considered 
as  constant  (see  Panofsky  and  Dutton,  1983).  In  this  layer,  both  mechanical  and  thermal  forcing 
affect  the  turbulence  and  the  variation  of  the  mean  wind  speed,  temperature,  and  humidity.  Monin 
and  Obukhov  [1954]  introduced  two  scaling  parameters  for  velocity  and  length  that  are  indep>endent 
of  height  in  the  surface  layer.  The  parameters  are  the  friction  velocity,  ,  and  a  length  L  that 
depends  only  on  the  fluxes  of  heat  and  moisture  and  the  friction  velocity. 

In  the  surface  layer,  the  vertical  gradient  of  a  conservative  property  i  is  related  to  its  vertical  flux 
Si  by 

—  = - — - <j)  ('>) 

p  ■  t  ■  u.(z  +  z„] 

where  p  is  the  atmospheric  density,  k  is  von  Karmen’s  constant  (0.4),  is  the  friction  velocity,  Zo 
is  a  surface  roughness  parameter,  and  O  is  a  stability  function.  Jeske  [1973]  proposed  that  <I>  take 
the  form  of  the  Monin-Obukhov  logarithmic-linear  model 

<I>  =  1  4-  (3) 

for  stable  conditions  (air  warmer  than  the  sea)  and  that  <I>  take  the  form  of  the  KEYPS  relation 
[Lumley  and  Panofsky,  1964] 

-  4/3^  =  1  (4) 

Li 


for  unstable  conditions  (air  cooler  that  the  sea).  In  these  expressions,  a  is  taken  to  be  5.2,  (3  is  taken 
to  be  4.5,  Zq  is  0.00015  meter,  and  I'  is  the  gradient  form  of  the  Morin-Obukhov  scaling  length 
corrected  for  stability.  Jeske  uses  an  empirical  profile  coefficient  T  to  relate  L'  to  physical  observ¬ 
ables  of  temperature  and  wind  speed: 


TaU^r 

giTa  -  Ts) 


(5) 


where  Ta  is  the  air  temperature  in  Kelvin,  is  the  sea  surface  temperatu.  ?  in  Kelvin,  g  is  the  accelera¬ 
tion  of  gravity  (9.8  m/s^),  and  U  is  the  wind  speed  in  m/s.  Taking  the  conservative  property  as  poten¬ 
tial  refractivity  ,  the  scaling  function  is 


N  = 

« 


k  •  ANp 


(6) 
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where  jyVp  is  the  potential  refractivity  difference  between  the  surface  and  a  reference  height,  Zi,  in 
meters.  The  universal  stability  function  xp  is  related  to  <1>  through 


= r  [ 


(7) 


The  vertical  modified  refractivity  profile  is  computed  as 

N. 


M  =  0.125  •  Z  + 


(8) 


The  surface  layer  model  of  Liu,  Katsaros,  and  Businger  [1979],  also  referred  to  as  the  LKB 
model,  is  based  on  simultaneously  solving  the  diabatic  profile  equations  for  velocity,  temperature, 
and  humidity,  given  by 

0  -0s  _  HZ/Zj)  -  xpQ  (9a) 

0,  1.14  •  k 


Q-Qs  ^  -  \pQ 

Q,  1.14  •  k 


(9b) 


U  -Us  _  ln(Z/Zo)  -  xpu 
U,  k 


where  0  is  the  potential  temperature  and  Q  is  the  specific  humidity. 
For  unstable  conditions,  the  stability  functions  are 


V'©  =  =  2  •  In 


1  +  (1  - 

2 


(9c) 


(10a) 


y’u 


=  2  .  +  ,„1M1-^A)-/L2  .  ,a„-,l-16Z/L)V4  +  f 


For  stable  conditions,  the  stability  functions  are 

V'0  =  V'q  =  V'l/  =  -1  ’ZjL 
The  Monin-Obukhov  length  is 


L  = 


(0  (1  +  0.61  •  0  •  vy- 


g  -  k  (0^  •  (1  +  0.61  •  0  +  0.61  •  0  •  e.) 

The  surface  roughness  parameter  Zq  is  related  to  the  wind  drag  coefficient  Qj  by 


Zq  =  10  '  exp 


,-i 


where  Q)  is  computed  for  the  wind  speed  at  10  meters,  which  is  given  by 


(10b) 

(10c) 

(11) 


(12) 
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(13) 


The  remaining  unknowns,  Zj  and  Zg,  are  related  to  the  Reynolds  roughness  number,  which  is 
related  to  the  surface  roughness  parameter  and  the  friction  velocity  [Liu,  Katsaros,  and  Businger, 
op.  cit.]. 


Equations  9  through  13  are  solved  iteratively  to  obtain  the  estimates  of  the  temperature,  humid¬ 
ity,  and  velocity  scaling  functions.  Equation  9  is  used  to  calculate  the  vertical  temperature  and 
humidity  profiles  and,  with  the  hydrostatic  approximation  to  allow  the  calculation  of  the  pressure 
profile,  equation  1  is  used  to  compute  the  vertical  profile  of  refractivity. 


The  LKB  model  is  in  a  sense  more  pleasing  than  the  Jeske  formulation  in  that  it  does  not  require 
an  empirical  connection  between! '  and  the  physical  observables.  The  added  complexities  of  solving 
three  equations  for  three  unknowns  are  generally  insignificant  except  in  highly  stable  conditions 
where  the  problem  is  ill-formed  and  the  iterations  fail.  There  are  probably  techniques  for  solving 
these  ill-formed  problems,  but  that  topic  is  beyond  the  scope  of  this  analysis. 


RESULTS 

METEOROLOGICAL  MEASUREMENTS 

The  values  of  Zq  and  L ',  derived  from  the  LKB  formalism,  are  listed  in  table  2.  Figures  2  through 
21  show  the  vertical  modified  refractivity  profiles  corresponding  to  the  entries  in  table  2.  Data  pro¬ 
cessed  by  the  PP-11  are  shown  as  open  circles;  data  processed  by  the  MRS  are  shown  as  crosses. 
The  solid  line  is  the  M  profile,  calculated  using  the  LKB  formalism,  for  the  surface  observations 
listed  in  table  2.  Estimates  of  the  modified  refractivity  vertical  gradient  {dMIdZ)  are  listed  in  the  last 
three  columns  of  table  2.  Gradients  derived  from  radiosonde  measurements  are  calculated  by  fitting 
a  linear  regression  line  to  the  PP-11  and  MRS  observations  in  the  height  interval  of  100  to  250 
meters.  The  gradients  derived  from  the  bulk  measurements  are  calculated  by  fitting  a  linear  regres¬ 
sion  line  to  the  LKB  M  profile  in  the  height  interval  of  50  to  100  meters.  Figure  22  compares  these 
M  gradients  showing  dMldZ  calculated  using  LKB  techniques  as  the  abscissa  and  dMIdZ  derived 
from  the  radiosondes  as  the  ordinate.  There  is  generally  good  agreement  between  the  gradients 
derived  from  upper-air  observations  and  the  gradients  computed  from  bulk  surface  measurements. 
There  are  six  cases  that  are  not  plotted  in  figure  22.  These  cases  are  marked  with  an  asterisk  in  the 
dMIdZ  columns  of  table  2.  On  the  days  of  16  and  17  September,  the  significant  discrepancies  in  the 
gradients  are  likely  attributable  to  errors  in  processing  by  the  PP-11  and  MRS.  Figures  4  (MRS  fail¬ 
ure)  and  7  (MRS  and  PP-11  failure)  clearly  illustrate  these  errors.  The  discrepancies  on  20  Septem¬ 
ber  (see  figure  11)  are  not  so  easily  explained,  but  may  have  resulted  from  inadequate  surface  mea¬ 
surements.  The  radiosonde  measurements  appear  to  show  a  surface  duct  about  10  meters  thick, 
whereas  the  surface  layer  model  shows  a  surface  duct  about  1  meter  thick.  However,  in  general,  there 
is  reasonable  agreement  between  dMIdZ  measured  by  radiosondes  and  dMIdZ  calculated  from  the 
surface  meteorological  measurements. 

Figures  23  through  42  compare  the  vertical  profiles  of  temperature,  relative  humidity,  and  modi¬ 
fied  refractivity  of  the  PP-11 -processed  radiosonde  data  with  the  vertical  profiles  calculated  from 
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the  bulk  surface  meteorological  measurements.  In  most  of  these  cases,  there  is  very  good  agreement 
between  the  gradients  of  the  vertical  profiles  (from  the  surface  to  a  height  of  75  meters) .  The  surface 
and  upper-air  temperature  profiles  around  0700  GMT  are  markedly  different  on  the  morning  of  22 
September,  as  shown  by  figures  37  and  38.  In  the  morning,  the  air  was  cool  (about  1 1 .6°  C)  but  rapid¬ 
ly  increased  in  temperature  (12.3°C  half  an  hour  later);  the  radiosonde  measurements  indicate  a  sur¬ 
face  inversion  layer.  Part  of  the  discrepancies  in  the  temperature  gradients  may  be  due  to  radiational 
heating  and  cooling  of  the  land  (wind  speed  was  low,  about  2  knots). 

Figure  41  illustrates  the  only  strongly  stable  surface  measurements  made  at  Gavres.  The  wind 
speed  was  low,  about  3  knots,  and  coming  from  the  land  (060  degrees).  It  is  likely  that  the  large  dif¬ 
ferences  in  the  temperature,  humidity,  and  M  profiles  are  due  to  radiational  heating  by  the  land. 

ANGLE-OF-ARRIVAL  ESTIMATION 

Figures  43  and  44  show  calculated  estimates  of  the  phase-front  angle-of-arrival  for  frequencies 
of  10.5  and  35  GHz.  It  is  assumed  that  the  transmitter  is  1 8 .6  meters  above  the  water,  that  the  receiver 
is  10  meters  above  the  water,  and  that  the  antennas  are  separated  by  9.8  km.  The  angle-of-arrival 
is  estimated  as  the  real  part  of  an  eigenvalue  solution  to  the  wave  equation  evaluated  at  the  surface 
of  the  ocean  (see  Budden,  1961).  In  essence,  the  eigenvalue  should  be  evaluated  at  the  receiver 
height,  but,  for  the  low  sited  antennas  considered  here,  an  evaluation  at  the  surface  should  be  an  ade¬ 
quate  approximation.  The  abscissa  is  the  propagation  factor  that  results  from  a  power  evaluation  of 
each  eigenvalue  or  mode.  It  includes  effects  of  the  height-gain  functions  and  is  a  measure  of  the 
power  contained  in  the  mode. 

Three  cases  are  shown  on  each  of  these  plots.  The  triangles  represent  estimates  of  the  an<ile-of- 
arrival  when  the  refractive  conditions  are  normal  (over  water  the  gradient  is  0. 13  M/m).  A  .normal 
atmosphere  is  considered  as  a  reference  atmosphere  in  propagation  analysis.  The  squares  represent 
angle-of-arrival  estimates  for  the  siuface-reff  active  environment  measured  at  1 1 30  GMT  on  2 1  Sep¬ 
tember,  but  with  an  additional  assumption  of  a  perfectly  smooth  ocean  surface.  The  circles  represent 
angle-of-arrival  estimates  for  the  surface-refractive  conditions  at  1 130  GMT  on  21  September,  with 
a  rough  ocean  surface  correspondirig  to  a  wind  speed  of  10  knots.  There  are  only  slight  differences 
between  rough  and  smooth  surface  conditions. 

The  significant  aspects  of  figures  43  and  44  are — 

1.  Propagation  assessments  using  the  measured  meteorological  conditions  imply  that  the 
angle-of-anrival  changes  from  about  1.5  mrad  (in  a  normal  atmosphere)  to  about  6.2 
mrad  for  the  specific  refractive  conditions  measured  at  1130  GMT  on  21  September. 

2.  The  changes  in  angle-of-arrival  are  approximately  independent  of  frequency,  which  is  in 
general  agreement  for  a  ray  optics  treatment  of  the  propagation. 

3.  Surface  roughness  has  little  or  no  impact  on  the  angle-of-arrival. 
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CONCLUSIONS 

TVo  weeks  of  surface  and  upper-air  meteorological  measurements  were  completed  in  Lorient, 
France,  in  support  of  an  international  effort  which  was  coordinated  by  NATO/AC243/Panel  3 
Research  Study  Groups  8  and  21.  These  propagation  and  meteorological  measurements  were 
designed  to  assess  the  effects  of  the  lower  atmosphere  on  the  propagation  of  the  electromagnetic 
phase  front  at  frequencies  from  10  to  94  GHz. 

The  numerical  modeling  of  propagation  from  waveguide  theory  at  10.5  and  35  GHz  in  normal 
and  ducting  environments  indicates  that  the  evaporation  duct  can  raise  the  angle-of-arrival  from 
about  1.5  mrad  to  about  6.2  mrad.  These  changes  in  angle-of-arrival  may  be  significant  for  some 
close-in  weapon  systems.  In  addition,  the  modeling  indicates  that  surface  roughness,  at  least  on  this 
short  path,  has  little  or  no  effect  on  angle-of-arrival. 

The  indications  of  angle-of-arrival  frequency  independence  should  be  examined  in  greater 
detail.  A  ray  optics  evaluation  and  comparison  to  the  waveguide  solutions  may  lead  to  additional 
insight  into  the  problem  of  propagation  assessment. 
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Figure  1.  The  geographical  region  around  Lorient,  France, 
in  the  Bay  of  Biscay. 
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Figure  2.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  15  September  1993, 0730  GMT. 
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15  September  1993,  1730  GMT 
L*  =  -82.02  tn 
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Figure  3.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  15  September  1993, 1730  GMT. 
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Figure  4.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  16  September  1993, 0715  GMT. 
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Figure  5.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  16  September  1993, 1100  GMT. 


16  September  1993,  1415  GMT 
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Figure  6.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  16  September  1993, 1415  GMT. 
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17  September  1993,  0715  GMT 
L*  a  -37.79  m 


300  310  320  330  340  350  360  370 

M-UNITS 


Figure  7.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  17  September  1993, 0715  GMT. 


1 7  September  1 993, 1115  GMT 
L*«-2.19in 


300  310  320  330  340  350  360  370 


M-UNITS 

Figure  8.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  17  September  1993, 1115  GMT. 
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17  September  1993, 1130  GMT 
L‘  =  -2.19m 


300  310  320  330  340  350  360  370 

M-UNITS 


Figure  9.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  17  September  1993,  1130  GMT. 


1 7  September  1 993, 1 41 5  GMT 
L*  =  -6.l7ni 


Figure  10.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  17  September  1993, 1415  GMT. 
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20  September  1 993, 1 200  GMT 
L*  =  -41.31  m 


Figure  11.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  20  September  1993, 1200  GMT. 


20  September  1993, 1430  GMT 
L*« -20.74  m 


Figure  12.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  20  September  1993, 1430  GMT. 
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21  September  1993,  0800  GMT 
L’  a  -44.65  m 


M-UNITS 


Figure  13.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  21  September  1993, 0800  GMT. 


21  September  1993, 1130  GMT 
L’a -155.6  m 


300  310  320  330  340  350  360  370 

M-UNITS 


Figure  14.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  21  September  1993, 1130  GMT. 
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21  September  1993, 1400  GMT 
L* » -72.39  m 


300  310  320  330  340  350  360  370 

M-UNITS 

Figure  15.  Vertical  proxies  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  21  September  1993, 1400  GMT. 


22  September  1993, 0700  GMT 
L’  ■  -0  m 


300  310  320  330  340  350  360  370 


M-UNITS 

Figure  16.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  22  September  1993, 0700  GMT. 
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22  September  1993,  0730  GMT 
L’s-l.OOm 


300  310  320  330  340  350  360  370 

M-UNITS 


Figure  17.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  22  September  1993, 0730  GMT, 


22  September  1993, 1100  GMT 
L’« -11.00  m 


300  310  320  330  340  350  360  370 

M-UNITS 


Figure  18.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  22  September  1993, 1100  GMT. 
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22  September  1993, 1400  GMT 
L’  =  -8.03  m 


Figure  19.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  22  September  1993, 1400  GMT. 

23  September  1993, 1130  GMT 
L’  =  +4.67  m 


Figure  20.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  23  September  1993, 1130  GMT. 
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23  September  1993, 1400  GMT 
L’ « -28.63  m 


Figure  21.  Vertical  profiles  of  the  modified  refractivity 
measured  by  upper-air  radiosondes  and  surface  layer 
calculations  for  23  September  1993, 1400  GMT. 


dM/DZ  FROM  SURFACE  OBSERVATIONS  (LKB) 


Figure  22.  Gradients  of  modified  refractivity 
computed  from  surface  layer  observations 
and  compared  to  gradient  measured  by  upper- 
air  radiosondes. 
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Location  Gaures,  France  (47  42  N,  03  20  U) 
Datc/Tine:  03-15-1333  07:18:25  GMT-2 


14  T  18  60  RH  100  320  M  350 


U(kt)=Z4  0  P(««b)=  930.9  T(C)=*16.1  RH=  84.0  Ts=»15.9  dP^dZ=-0.117 

Ustar=«6.423E-01  Tstar=*7.079E-e3  Qstar=-7.924E-05  20=*4 . 552E-04  L' =-4.447E«03 
Zstap=*1.127E*01  Delta=*9.475E*0O 

Figure  23.  Vertical  profiles  of  temperature,  relative  humidity,  and  modified 
refractivity  measured  by  radiosondes  (PP-11  processor)  compared  to  profiles 
computed  from  surface  layer  observations  for  15  September  1993, 0730  GMT. 

Location  :  Gawes,  France  (47  42  N>  03  20  U) 

Date/Il«ic:  09-15-1393  17:27:03  (31T-Z 


U(kt)=10.O  P(iib)=100i.9  T(C)=*lb.3  RH=  62.0  Ts=*16.6  <lP/d2=-0.ll8 

Ustar=*2.357E-01  Tstar=-1 .309E-0Z  Qatar =-2.098E-O4  Z0=*1.535E-®4  L’ =-B.20ZE*01 
Zstar=*1.507E*01  Delta=»1.313E»81 


Figure  24.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  15  September  1993, 1730  GMT. 
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Location  ;  Gaures,  France  (17  12  N.  03  20  U) 
Date/Tine;  09-16-1993  07:1O;3O  Gf1T-2 


U(kt)=  8.0  P(nb)=1009.9  T(C)=*11.7  RH=  71.0  Ts=*15.5  dP/dZ=-O.120 

Ustar=*1.811E-01  Tstar=-3.812E-02  Qstar=-1.501E'01  20=*l.e31E-01  L' =-3.723E*01 
2star=*9.312E*0e  Delta=*8.220E*e0 


Figure  25.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  16  September  1993, 0715  GMT. 


Location  :  Gaures,  France  (17  12  N,  03  20  U) 
Date/Tine:  09-16-1993  10:59:26  KlT-2 


U()ct)=  8.0  r(nb}=1011.9  T(C)=*15.7  »i=  76.0  Ts=*16.7  dP/'dZ=-0 . 120 

Ustar=*1.816E-01  Tstar=-1.819E-02  Qstar=-1.337E-01  20=»l. 0352-01  L’ =-3.381E«01 
Zstar=*a.228E*0e  DeIta=*7.Z10E«e0 


Figure  26.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  16  September  1993, 1100  GMT. 
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Location  ;  Gaures.  France  (47  42  11,  03  20  U) 
Oate/TiMc;  03-16-1393  14:17:26  GnT-2 


U(kt)=  9.0  P(iib)=1013.9  T(C)=*15.9  M1=  71.0  Is=*lb.7  dP/dZ=-0.1Z0 

Ustar=>2.082E-61  Tstar=-3.6BSE-eZ  Qstar=-1.443E-04  Z0=*1 .266E-04  L'=-5.123E*OI 
Zstar=«8.913E*Oe  Delta=*8.660E*e0 


Figure  27.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  16  September  1993, 1415  GMT. 


Location  :  Gavres,  France  (4?  42  N,  03  ZO  U) 
Date/time:  09-17-1993  07:15:52  GHT-Z 


U(kt)=  Z.O  P(Mb)*1017.9  T(C)=  *9.0  IIH=  79.0  Ts=*15.1  dP/'d2=-0.123 

Ustar=*4.797E-0Z  tstar=-3.875E-01  qstar=-Z.754E-04  ZO=*3.550E-O5  L' =-3.779E-01 
Zstar=*2.75ZE*0e  Delta=*2.545E*O0 


Figure  28.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  17  September  1993, 0715  GMT. 
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Location  :  Gaurcs,  France  (17  12  H,  03  20  U) 
Date/'Tine:  09-17-1393  11:16:07  GtlT-Z 


U(kt)=  3.0  P(nb)=:iei7.9  T(C)=*15.8  »H=  57.0  ts=*17.5  dP/dZ=-0.120 

Ustap=*6.31BE-0Z  Tstar=-8.753E-0Z  Qstar=-2.b91E-01  Z0=*2.b62E-G5  L' =-2.1B9E»G0 
Zstar=*5.68BE*00  Delta=*5.000E«00 


Figure  29.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  17  September  1993, 1115  GMT. 


Location  :  Gavres,  France  (17  12  It,  03  20  U) 
Date/Tine:  09-17-1993  11:39:38  GOT-Z 


U(kt)=  3.0  P(nb)=1017.9  T(C)=*15.8  RH=  57.0  Ts=*17.5  dP/dZ=-0.1Ze 

Ustar=*6.318E-02  rstar=-B.753E-0Z  Qstar=-Z.691E-©1  Z0=«2.662E-05  L' =-2.189E«00 
Zstar=«5.b88E'K)e  Delta=*5.080E*00 


Figure  30.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  17  September  1993, 1130  GMT. 
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Location  :  Gaures.  France  (4?  tZ  N<  03  20  U) 
Date/Tine:  0S-17-1S93  M:1Z:1Z  GflT-Z 


U(kt)=  5.0  P(nb)=ieH.9  T(C)=*16.2  RH=  62.0  Ts=*lB.O  dP/dZ=-0 . 120 
Ustar=*1.06«-01  Tstar=-9.069E-02  Qstar=-2.507E-Od  20=*4 . 18BE-05  L' =-6 . 16BE*O0 
Zstar=*7.538E*00  0elta=*6.63SE*e0 


Figure  31.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  17  September  1993, 1415  GMT. 


Location  ;  Gaures.  France  (47  42  N,  03  20  U) 
Date/Tlne:  09-20-1993  12:14:51  GnT-2 


U(kt)=  5.0  P(Mb)=101S.9  T(C)=*16.6  l)H=  94.0  Ts=*17.0  dP/dZ=-0.120 

Ustar=*l.O30E-Ol  Istar=-l.BS8E-02  Qstar=-l.7B2E-06  20=*4 . 188E-05  L' =-4.131E«01 
Zstar=-2.553E-01  Delta=-9.990E*01 


Figure  32.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  20  September  1993, 1200  GMT. 
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Location  Gawres,  France  (47  42  H,  03  20  U) 
Date/Tine:  03-20-1993  14:17:41  GnT-2 


U(kt)=  5.0  P(nb)=iei5.9  T(C)=«16.4  BH=  66.0  Ts=*17.0  dP/dZ=-e . 120 

Ustar=*l.e38E-01  Tstar=-2 .aSlE-02  Qstar=-5.521E-05  20=*4 . 18BE-05  L' =-2.074E»01 
Zstar=*3.714E*0e  Delta=«3.ZZ5E*00 


Figure  33.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  20  September  1993, 1430  GMT. 


Location  :  Gawr«s,  France  (47  42  N,  03  20  M) 
Date/Tine:  09-21-1993  08:02:12  GUT-E 


H(kt)=  7.0  P(nbl=1014.9  T(C)=*15.9  »i=  82.0  Ts=»16.4  dP/'d2=-0.120 

Ustar=*1.533E-01  Tstar=-2.468E-02  Qstar=-7.777E-05  Z0=*7.53BE-O5  L’ =-4.465E*01 
2star='»b.019E'»00  Delta=*5.230E*00 


Figure  34.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  21  September  1993, 0800  GMT. 
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Location  :  Gaures,  France  (■i?  12  N,  03  20  U) 
Date/'T Ine:  03-21-1333  11:25:12  011-2 


14  T  18  60  RH  100  320  M  360 


U(kt)=10.0  P(iib)=1017.3  7(C)=*17.0  RH=  71.0  Ts=«l7.2  dP/d2=-e .  120 

Ustar=«2.350E-Ol  Tstar=-a.Sl2E-03  Qstar=-l.e05E-01  Z0=«1 .52BE-01  L' =-1.556E»02 
2star=«1.015E*01  Deita=*8.675E*00 

Figure  35.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  21  September  1993,  1130  GM. 

Location  :  Gawres,  France  <47  12  N,  03  20  U) 

Date/Tine:  03-21-1333  14:06:17  GflT-E 


14  T  18  50  RH  100  320  M  360 


U(kt)=  3.0  P(nb)=1017.3  I(O=*17.0  »»=  73.0  Ts=*17.5  dP/d2=-e.l20 

Ustar=*2.076E-01  tstar=-2.276E-92  Qstar=-1.202E-04  Z0=*1.261E-04  L' =-7.233E*01 
Zstar=*3.613E*0e  Delta=«6. 3002*00 

Figure  36.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface  layer 
observations  for  21  September  1993, 1400  GMT. 
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Location  :  Gawres,  France  (47  42  M,  03  20  U) 
Date/Tioe:  03-22-1^83  07:04:30  GnT-2 


U(kt)=  Z.O  P(i»h)=1012.9  T(C)=»11.6  RH=  84.0  Ts=*15.6  dP/dZ=-©.121 

Ustar=*4.6aOE-OZ  Tstar=-2.436E-01  Qstar=-Z.051E-04  20=*3.b3aE-05  L'=-5.667E-01 
Zstar=*Z.667E*0e  Delta=*Z.4Z0E*e0 


Figure  37.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  22  September  1993, 0700  GMT. 


Location  :  Gaores,  France  (47  42  N,  03  20  U) 
Date/'TiMe;  03-22-1993  07:28:20  GFIT-Z 


U(kt)=  2.5  P(i»b)=1811.9  T(C)=»12.3  BH=  83.0  Ts=*X5.7  dP^dZ=-8.121 

Ustar=*5.544E-02  Istar=-1.909E-81  Qstar=-1.807E-G4  Z0=*3.043E-O5  L' =-l .803E*00 
Zstar=*3.ei6E*00  Delta=*Z.730E*00 


Figure  38.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  22  September  1993, 0730  GMT. 
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Location  :  Gaures,  France  (47  4Z  N,  03  ZO  U) 
Date/Tine:  09-ZZ-1933  11:15:Z3  GfIT-Z 


U(kt)=  b.e  P(iib)=iei5.9  I(C)=«15.5  RH=  65.9  Ts=«17.0  dP.^dZ=-© .  1Z9 

Ustar=»1.300E-01  Tstar=-7.577E-0Z  Qstap=-Z .e77E-01  20=«5.7Z7E-05  L' =-l . 100E*01 
2star=*7.977E*0e  Delta=*7.035E*e0 


Figure  39.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  22  September  1993, 1100  GMT. 


Location  :  Gavres,  France  (47  4Z  N,  03  ZO  U) 
Date/Tine:  09-Z2-1993  14:08:57  GHI-Z 


UOct)=  5.0  P(nb)=iei5.4  T(C)=*16.e  60.0  Ts=*lB.0  dP^dZ=-e.lZ0 

Ustar=*l.eS6E-01  Tstar=-5.950E-OZ  Qstap=-Z.436E-04  ZO=*4.188E-05  L‘ =-B.O30E*OO 
Zstar=«8.177E*00  Delta=«7.160E*00 


Figure  40.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  22  September  1993, 1400  GMT. 
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Location  :  Gawres,  France  (47  42  N,  03  20  U) 
Date/Tine:  03-23-1993  11:17:11  GHI-2 


U(Jft)=  3.0  P(*ib)=iei6.9  I(C)=*18.5  RH=  56.0  Is=*16.B  dPxdZ=-0. 119 

Ustar=*5.216E-02  Tstap=*6.ae3E-02  Qst*r=-1.427E-©4  20=*3.208E-05  L' =*4 ,676E*0G 
Zstap=-9.99eE*01  Delta=-9.99OE»01 


Figure  41.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  23  September  1993, 1130  GMT. 


Location  :  Gawres.  France  (47  42  N,  03  20  U) 
Date/Tine:  09-23-1993  14:08:48  aiT-2 


U(kt)=  8.0  P(nb)=iei6.4  T(C)=*16.B  RH=  73.0  Ts=*18.0  dP/dZ=-e.lZ0 

Ustar=*1.8ZeE-01  Tst*r=-5.0O1E-O2  Qstar=-1.546E-04  ZO=»1.03aE-O4  L' =-Z.B63E*01 
Zstar=*8.696£«0e  I>elta=*7.5B0E*e0 


Figure  42.  Vertical  profiles  of  the  modified  refractivity  measured  by 
upper-air  radiosondes  compared  to  profiles  computed  from  surface 
layer  observations  for  23  September  1993, 1400  GMT. 
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Figure  43.  Estimated  angle -of-arrival  at  10.5  GHz  for 
1130  GMT  on  21  September  1993. 


RECEIVED  POWER  CONTRIBUTION  (dB) 


Figure  44.  Estimated  angle-of-arrival  at  35  GHz  for 
1130  GMT  on  21  September  1993. 
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